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SUMMARY 


V 

An  -investigation  Is  carried  out  into  the  structure  of  the  laminar  boundary 
layer  originating  from  the  forward  stagnation  point  of  a  prolate  spheroid  at 
Incidence  in  a  uniform  stream,  assuming  that  the  external  velocity  distribution 
Is  given  by  attached  potential  theory.  The  principal  new  results  of  the 
study  are  ~ 

(()  A  new  transformation  of  the  body  coordinates  is  devised  which 

facilitates  the  computation  of  the  solution  near  the  nose^ 

(  '  ' 4 1 if-  Two  variations  of  the  standard  box  method  of  solving  the 

equations  are  devised  to  enable  solutions  to  be  computed  in  regions  of 

cross-flow  reversal.  They  are  referred  to  as  the  zig-zag  box  and  the 
y  \ 

characteristic  box.  J 

4*444-  Whereas  in  two  dimensional  flows  the  effect  of  the  boundary 
layer  approaching  separation  on  the  external  flow  may  be  represented  by 
a  blowing  velocity.  In  the  present  study  we  find  that  this  is  only  true 
near  the  windward  line  of  symmetry^  Near  the  leeward  line  of  symmetry 
the  blowing  velocity  must  be  replaced  by  a  suction  velocity  even  though 
the  boundary  layer  is  being  significantly  thickened. 

The  boundary  layer  over  the  whole  of  the  spheroid  cannot  be  computed 

in  an  integration  from  the  forward  stagnation  point.  If  o  <  6°,  the  access- 

„  _ _ •* 

ible  region  is  largely  bounded  by  the  separation  line  and  develops  a  wedge-like 


shape  whose  apex  is  named  the  ok  of  accessibility,  pointing  towards  the 
nose  of  the  spheroid.  On  the  windward  side  of  this  line  the  solution 
develops  a  singularity;  on  the  leeward  side  the  situation  Is  less  clear 
but  It  is  believed  to  occur  there  too. 
x  ^  4^  For  a  >  15^  the  accessible  region  on  the  leeward  side  of 
the  ok  is  largely  determined  by  the  external  streamline  through  the  ok. 
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I.  INTRODUCTION 

The  determination  of  the  three-dimensional  boundary  layer  on  a  body  when 
subject  to  a  prescribed  pressure  gradient  has  attracted  much  Interest  during 
the  last  twenty  years  since  the  general  principles  which  must  govern  any 
sound  numerical  approach  were  laid  down  by  Raetz  (1957).  These  follow  from 
the  realization  that  the  momentum  equations  are  diffusive  In  the  direction 
normal  to  the  body  and  wave- like  In  planes  parallel  to  the  body,  the  direction 
of  propagation  being  along  the  local  stream  direction.  Since  in  general  this 
direction  varies  across  the  boundary  layer  it  is  possible  to  identify  zones 
of  Influence  and  of  dependence  for  any  point  P  of  the  boundary  layer  each 
being  a  curvilinear  wedge  with  generators  composed  of  normals  to  the  body 
and  having  the  normal  through  P  as  vertex  line.  The  computation  of  the 
boundary  layer  properties  at  P  must.  In  some  sense,  make  use  of  Its  properties 
at  all  points  In  the  zone  of  influence  on  P.  The  same  remark  applies  to  two 
dimensional  boundary  layers,  the  zone  of  dependence  now  reduces  to  a  plane. 

The  complication  in  three-dimensional  flow  Is  that  this  zone  can  be  very 
broad,  especially  near  separation  where  the  skin-friction  direction  may  be 
inclined  at  an  obtuse  angle  to  the  main  stream  and  even  worse  situations  may 
occur.  Beyond  separation  the  zone  of  dependence  in  two  dimensions  is  a  plane 
extending  both  upstream  and  downstream  In  terms  of  the  mainstream  while  In 
three-dimensions.  It  might  extend  in  all  directions,  at  least  near  P. 

The  concept  of  separation  and  Its  relation  to  that  of  accessibility  is  of 
some  subtlety  and  requires  careful  consideration.  The  problems  which  can  arise 
have  been  extensively  studied,  notably  by  Maskell  (1955),  Llghthill  (1963)  and  Wang 
(1976), the  last  of  whom  also  gives  a  well-balanced  review  of  previous  work. 
Accessibility  Is  defined  as  follows.  A  point  P  of  the  boundary  layer  Is 
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said  to  be  accessible  from  the  forward  stagnation  point  0  If  the  velocity 
field  at  P  can  be  computed  In  terms  of  the  Initial  conditions  at  0  and 
the  boundary  conditions  on  the  body  and  In  the  external  stream.  The  boundary 
of  this  region  consists  of  the  normals  to  the  body  through  a  curve  on 
the  body,  t ^  being  often  closed  and  encircling  0  but  It  could  extend  to 
Infinity.  Separation  Is  a  line  drawn  on  the  body  which  forms  the 
boundary  of  all  limiting  streamlines  emanating  from  0.  If  Is  closed, the 
normals  to  It  bound  the  region  of  accessibility  but  to  Identify  and 
requires  additional  arguments,  although  Wang  (1976)  noted  a  general  agreement 
that  this  Is  always  the  case.  Later  on  however  we  shall  demonstrate  that 
Important  exceptions  can  occur. 

It  has  been  argued  that  ts  Is  a  limiting  streamline,  l.e.  a  skin-friction 
llne.passlng  through  Isolated  singular  points  of  the  solution  (Llghthlll,  1963) 
or  Is  an  envelope  of  limiting  streamlines  (Maskell,  1955).  There  is  a  signifi¬ 
cant  difference  between  the  nature  of  the  boundary- layer  solution  according  to 
which  option  occurs.  For  if  it  Is  an  envelope, the  skin  friction  component  in 
a  direction  perpendicular  to  £$  must  have  an  algebraic  singularity  at  £$ 
with  index  probably  equal  to  1/2.  Further,  the  displacement  thickness  Is  also 
singular  at  £$  and  so  the  solution  terminates  at  i$.  On  the  other  hand  If 
£$  Is  a  skin-friction  line  then  It  Is  possible  for  the  solution  to  be  smooth 
at  £$  and  so  may  be  continued  beyond  the  region  of  accessibility  from  the 
forward  stagnation  point  although,  to  be  sure,  additional  information  must 
then  be  supplied  to  specify  It  uniquely.  Now  in  real  flows  ore  can  reasonably 
expect  that  the  pressure  gradient  adjusts  near  to  prevent  the  singularity  - 
otherwise  a  contradiction  might  well  occur  since  the  corresponding  singularity 
In  the  displacement  thickness  provokes  a  singularity  in  the  pressure  gradient 
forcing.  In  turn,  separation  to  occur  earlier.  This  argument  has  already  been 


Identified  this  with  the  free  vortex-layer  type  of  separation  previously 
described  by  Maskell  (1955).  Another  possible  type  of  separation  on  smooth 
bodies,  first  Identified  by  Howarth  (1951),  Is  a  collision  between  two 
boundary  layers  originating  from  different  parts  of  the  flow  field, but  for  a 
long  time  such  a  phenomenon  was  regarded  as  exceptional.  Recently,  however, 
a  number  of  additional  studies  have  been  reported  (see  Stewartson,  Cebecl 
&  Chang  1980,  for  examples)  and  it  seems  worthwhile  to  bear  It  in  mind  as  a 
possible  termination  of  the  calculations.  The  novel  feature  is  that  the 
forward  calculation  of  the  one  boundary  layer  develops  unexceptional ly 
until  It  collides  with  the  other,  the  exact  position  of  the  collision  being 
determined  by  mechanisms  which  are  not  yet  clear  but  which  probably  depend 
on  balancing  the  forward  momentum  In  each  of  the  boundary  layers. 

This  discussion  indicates  that  understanding  of  these  boundary  layers, 
even  when  the  pressure  gradient  Is  prescribed, is  still  Incomplete.  Our  aim 
in  this  report  is  partly  to  offer  an  alternative  numerical  method  which  will 
hopefully  enable  us  to  carry  the  Integration  further  towards  and 
clarify  some  of  the  outstanding  issues.  In  addition  we  have  some  new  goals. 

It  is  well-known  that  in  problems  of  engineering  Interest  the  majority 
of  the  boundary  layer  Is  turbulent,  transition  occurring  near  the  nose  of 
the  body  at  high  angles  of  Incidence.  In  order  to  study  these  flows  we 
first  need  an  accurate  method  of  computing  the  laminar  flow  In  this 
neighborhood.  The  method  used  by  Wang  (1976)  has  limitations  here  because 
he  does  not  choose  to  remove  the  singularities  In  the  metric  parameters  at 
the  nose  and  while  the  use  of  coordinates  based  on  properties  of  the  Imposed 
mainstream  formally  removes  them  (Gelssler  1975),  they  lead  to  additional 
difficulties  [Cebecl,  Kaups,  Moslnskls  &  Rehn  1973].  Recently  Cebecl, 
Khattab  &  Stewartson  (1980)  have  devised  a  set  of  body  coordinates  with 


extensively  used  In  two-dimensional  flows, culminating  In  the  Sychev-Smlth 
theory  of  separation  on  a  circular  cylinder  and  other  smooth  bodies  (Sychev 
1967,  1972,  Smith  1977a,  1979). In  which  the  flow  field  Is  regular  as  the 
skin  friction  changes  sign.  It  Is  noted  that  some  studies  of  the  analogous 
triple-deck  theory  In  three-dimensional  flows  have  already  been  made 
(Burggraf  1978,  Smith,  Sykes  and  Brighton  1977,  Sykes  1979)  and  the  flow 
near  separation  Is  again  regular.  Now  there  are  many  photographs  available 
of  three-dimensional  separation;  Wang  (1976)  has  included  a  selection  in  his 
review,  for  example,  and  Han  and  Patel  (1979)  have  published  some  for  laminar 
flows  past  spheroids.  These  strongly  suggest  that  &s  is  an  envelope,  but  for  the 
above  reason  we  are  reluctant  to  accept  that  this  is  really  the  case  and  are 
of  the  opinion  that  a  more  detailed  study  would  show  it  to  be  a  skin-friction 
line  as  Lighthill  suggested. 

There  have  been  a  number  of  careful  calculations  of  three-dimensional 
boundary  layers  in  the  past  particularly  by  Wang  (1970,  1974a, b,c,  1972,  1975) 
and  by  Cebeci  and  his  associates  (1973,  1976), but  for  all  the  computation 
of  the  separation  line  has  proved  difficult.  If  not  impossible.  The 
most  successful  are  due  to  Wang  for  boundary  layers  on  a  prolate  spheroid 
at  incidence;but  even  he  was  not  able  to  comment  usefully  on  the  nature  of 
the  solution  near  separation, and  the  computation  on  the  leeward  side  of  the 
body  proved  very  difficult.  He  took  careful  note  of  the  zones  of  dependence 
in  formulating  his  numerical  algorithm  and  used  an  implicit  finite-difference 
method  based  on  the  Crank-Nicholson  scheme.  For  an  angle  of  attack  a  =  30°  he 
concluded  that  the  closed  separation  line  i$,  which  we  have  already  discussed 
at  length,  is  replaced  by  an  open  separation  line  iQ,  which  starts  about  10% 
of  the  way  to  the  back  of  the  body  and  about  140°  from  the  windward  line  of 


symmetry,  the  limiting  streamlines  approaching  it  from  both  sides.  He 


respect  to  which  the  governing  equations  are  regular  everywhere  upstream  of 

Here  we  shall  use  this  system  of  coordinates  and.  Indeed,  make  a 
further  transformation  to  convert  the  formal  computation  of  the  flow  near 
the  nose  into  a  straightforward  numerical  problem,  so  that  we  shall  later 
be  In  a  strong  position  to  examine  the  stability  of  the  laminar  boundary  layer. 

Second  we  wish  to  have  a  method  easily  adaptable  to  the  Integration  of 
turbulent  boundary  layers.  For  such  flows  It  Is  known  that  a  variable  mesh 
normal  to  the  body  Is  desirable,  to  cope  with  this  characteristic  double- 
structured  form.  A  generalization  of  the  Cebeci -Smith  eddy-viscosity  model 
(Cebeci  &  Bradshaw  1977)  has  quite  good  prospects  for  providing  an  adequate 
description  of  such  flows  and  the  corresponding  equations  have  close 
similarities  to  those  for  laminar  flow.  This  requirement  can  be  met  by  using 
the  Keller-box  method  (Cebeci  &  Bradshaw  1977)  in  which  the  governing  equa¬ 
tions  are  reduced  to  five  first-order  finite-difference  equations  in  terms  of 
values  either  at  the  corners  or  the  centers  of  the  sides  of  the  box  to  give  a 
second-order  accurate  solution  at  the  center  of  the  box.  These  boxes  are  now 
stacked  one  on  top  of  the  other  with  one  set  of  edges  along  the  normals  to  the 
surface.  The  nonlinear  difference  equations  are  linearized  by  using  Newton’s 
method  and  then  solved  by  a  standard  procedure  to  give  the  values  of  the  five 
dependent  variables  on  the  fourth  normal  in  terms  of  the  other  three.  In  this 
way  the  solution  can  be  advanced  from  either  the  windward  or  leeward  line  of 
symmetry,  where  it  can  be  found  independently  downstream  from  the  nose.  We 
shall  refer  to  this  method  as  the  standard  box. 

A  deficiency  in  the  standard  box  is  that  it  ignores  Raetz's  principle  of 
the  zone  of  dependence  and  so  the  user  is  not  surprised  that  it  eventually 
breaks  down  as  separation  is  approached  and  particularly  after  the  circumferential 
skin-friction  has  changed  sign.  Two  modifications  are  then  brought  into  use. 
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First  Krause  et  al.  (1968)  suggested  that  a  zig-zag  scheme  might  be  helpful  in 
this  situation:  in  an  integration  advancing  from  the  windward  line  of  symmetry 
this  amounts  to  replacing  one  of  the  edges  of  the  box  on  the  upstream  (and 
known)  side  by  the  next  edge  further  to  leeward.  We  investigate  its  usefulness 
here,  referring  to  the  modification  as  the  zig-zag  box,  first  used  by  Cebeci 
(1979)  for  unsteady  flows. 

However  this  approach  does  not  entirely  mimic  the  wave-like  character  of 
the  momentum  equations.  For  this  the  box  should  be  replaced  more  appropriately 
by  a  curvilinear  surface  consisting  of  the  (projections  of  the)  streamlines  on 
planes  parallel  to  the  surface  passing  through  the  normal  edge  along  which  we 
wish  to  compute  the  values  of  the  dependent  variables.  For  this  purpose  we 
use  a  new  numerical  algorithm  developed  by  Cebeci  and  Stewartson  (1977),  which 
incorporates  this  requirement:  we  shall  refer  to  it  subsequently  as  the 
characteristic  box,  thus  acknowledging  that  the  streamlines  are  characteristics 
of  the  momentum  equations. 

With  these  schemes  at  our  disposal  we  integrate  the  boundary- layer  equa¬ 
tions  for  a  prolate  spheroid  for  a  thickness  ratio  t  =  1/4,  set  at  an  angle 
of  incidence  a  to  an  oncoming  stream  where  a  *  3°,  6°,  15°,  30°,  taking 
the  imposed  pressure  gradient  to  be  that  given  by  classical  potential  theory. 

In  parenthesis  we  note  that  the  experimental  pressure  variations  agree  closely 
with  this  theory  except,  significantly,  near  separation  (Patel  &  Choi  1979), 
(Meier  &  Kreplln  1979).  We  find  no  difficulty  in  determining  the  solution 
near  the  nose  and  up  to  the  line  of  cross-skin  friction  reversal  using  the 
standard  box.  Further  downstream,  both  the  zig-zag  box  and  the  characteristic 
box  were  used  and  both  found  to  be  effective,  particularly  on  the  leeward 
side,  where  there  is  an  extended  region  before  separation  occurs;  it  appears 
that  the  characteristic  box  is  the  more  accurate  and  so  in  the  final  presenta¬ 
tion  of  results  it  is  preferred. 
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The  most  complete  of  the  earlier  studies  of  the  flow  fields  Is  for 
a  =  6°,  due  to  Wang  (1975),  who  concludes  that  It  may  be  computed  as  far  as 
the  separation  line  £$,  which  has  the  shape  of  a  tongue  half-way  between 
the  windward  line  of  symmetry  i  and  the  leeward  i  and  extends  significantly 
towards  the  nose.  Our  results  indicate  the  same  general  trends,  but  the  tongue 
is  much  more  pronounced  and  takes  on  a  wedge-like  character,  which  we  shall  refer 
to  as  an  ok  (lit.  arrow,  Turk.).  The  boundary  of  the  region  of  accessibility, 
defined  by  the  curve  thus  largely  coincides  with  when  a  =  6°.  In 
addition  we  are  able  to  carry  the  computation  near  enough  to  this  line  on 
the  windward  side  to  provide  convincing  evidence  that  the  solution  develops 
a  singularity  there,  probably  of  the  type  described  by  Brown  (1965).  In 
contrast  to  Wang's  results  we  find  that  the  circumferential  skin-friction 
does  not  vanish  on  £$,  and  nor  does  the  tangential  component  along  On 

the  leeward  side  we  find  marked  quantitative  differences  in  the  magnitude  of 
the  skin  friction  near  £$  and  indeed  in  the  position  of  However,  the 

leeward  portion  of  £s  cannot  be  fixed  so  definitely  as  the  windward  side 
since  we  are  unable  to  reproduce  the  Brown  singularity  at  all;  nevertheless 
we  are  fairly  confident  that  the  position  of  £$  inferred  from  our  computations 
is  close  to  being  correct. 

At  a  =  15°  the  identification  between  i ^  and  £$  is  only  possible 
on  the  windward  side.  The  form  of  £$  is  similar  to  that  for  a  =  6° 
except  that  it  originates  at  s  4  -0.4  and  in  particular  the  solution  is 
clearly  singular  there.  On  the  leeslde, however,  i ^  is  determined  by  the 
external  streamline  passing  above  the  most  forward  point  of  the  windward 
separation  line.  The  reason  is  that  downstream  of  this  line  the  solution 
in  the  outer  part  of  the  boundary  layer  is  historically  dependent  on  flow 
properties  along  streamlines  which  have  passed  over  the  windward  line  of 


separation  and  these  properties  are  not  fully  known.  Thus  an  ok  forms  in  also 
when  a  =  15°,  but  for  a  different  reason,  and  now  there  is  no  possibility  of  an 
accurate  determination  of  the  leeside  shape  of  with  the  boundary  values  and 
initial  conditions  supplied.  In  fact  the  concept  of  separation  is  strictly 
meaningless  here  and  the  notion  of  open  separation  introduced  by  Wang  (1975) 
is  irrelevant,  although  no  doubt  important  in  many  practical  flows.  Similar 
remarks  apply  to  the  solution  when  a  =  30°. 

An  important  aspect  of  the  boundary  layer  is  the  role  it  plays  in  correcting 
the  external  velocity.  In  two-dimensional  studies,  this  may  be  interpreted 
as  an  injection  velocity  from  the  body  into  the  inviscid  stream,  especially 
near  separation.  The  corresponding  result  for  the  three-dimensional  boundary 
layers  under  investigation  here  is  the  same  over  the  majority  of  the  flow 
field  computed, but  near  i  beyond  the  reversal  of  circumferential  velocity 
the  injection  velocity  is  replaced  by  a  suction  velocity.  Although  this 
unexpected  result  may  easily  be  inferred  from  previous  studies,  it  does 
not  seem  to  have  been  noted.  There  is  a  tendency  for  the  boundary  layer  to 
thicken  in  this  region  but  nevertheless  its  effect  on  this  equivalent  injec¬ 
tion  velocity  is  the  opposite  to  that  expected. 
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II.  FORMULATION 

The  studies  in  this  paper  extend  those  of  Cebeci  et  al.  (1980)  and  we 
shall  largely  follow  their  notation  and  formulation.  The  governing  boundary- 
layer  equations  for  an  incompressible  laminar  flow  in  a  curvilinear  orthogonal 
coordinate  system  appropriate  to  a  prolate  spheroid  at  incidence  (Fig.  1)  are 


Figure  1.  Notation  for  prolate  spheroid  at  incidence. 
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where  (u,w)  are  the  velocity  components  parallel  to  the  body  surface  and 
in, respectively, the  meridional  and  azimuthal  directions,  while  v  is  the 


component  normal  to  this  surface.  Further  h^ ,  ah2  are  metric  coefficients, 
where  n^  and  h2  are  defined  by 


1  ♦  g2(t2-l) 


i  -  r 


1/2 


,  =  t(l  -  r) 


2vl/2 


(2.4) 


t  denotes  the  thickness  ratio  (=b/a)  of  the  elliptic  profile  and  £  =  x/a. 
The  parameter  K2/a  is  the  geodesic  curvature  of  the  surface  lines  £  =  const, 
with  ^  9iven  by 


h1h2(l  -  £2)1/2 


(2.5) 


The  solution  of  the  system  (2.1)  to  (2.5)  requires  boundary  conditions  and 
initial  conditions.  The  boundary  conditions  are: 

y-0,  u  =  v  =  w  =  0, 

.  ,  ,  %  (2.6) 

y  u  ue(x,e),  w  -*•  we(x,e) 

The  velocity  components  ug  and  wg  can  be  obtained  from  inviscid  theory, 
(Hirsh  and  Cebeci  1977),  being  given  by 

ue 

u“  *  V0(t)  cosa  cose  -  Vgo(t)  sina  sine  cose  (2.7a) 

00 

W 

-u-  s  VgQ(t)  sina  sine.  (2.7b) 


Here  e  denotes  the  angle  between  the  line  tangent  to  the  elliptic  profile 
and  the  positive  £-axis:  it  is  given  by 
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V'.'fr 


cose  = 


[i  +  As  - 1)],/2 


The  parameters  VQ(t)  and  VgQ(t)  are  functions  of  t,  defined  by 


Vt>  = 


(1  -  t‘) 

z — 


2,3/2 


t2)  -f-  In 


Trrrm 

[i 


n 


Vg0(t) 


2V0(t) 
2V0(t)  -  1 


As  formulated,  the  governing  equations  are  singular  at  the  nose 
(£  *  -1)  and  this  may  be  removed,  as  explained  in  Cebecl  et  al  (1980),  by 
the  following  transformation.  We  define  new  velocity  components  U,  W,  V 
by 

u  =  U  cose  +  W  sine,  w  =  W  cose  -  U  sine,  tv  =  V 
and  new  coordinates  X,  Y,  Z  by 


where 


X  =  s  cose,  Z  =  s  sine,  tY  =  y 


dS 

5“  = 


[!  +  c2(t2  -  1)1 


t(l  -  cF) 


1/2 


Then  (2.1)  to  (2.3)  reduce  to 


N  (lr+  ly)+  TV  ~L(UX  +  WZ)  =  0  . 


n(u  Iy  +  w  |y)  +  LW(WX  -  UZ)  +  V  !y  =  +  V  ^  ,  (2.14) 

3  Y 

N(u  |y  +  w  |y)  -  LU(WX  -  UZ)  +  V  |y  =  e2  +  v  ^2  *  (2,15) 

^  3  Y 

where ^  and  ^  are  pressure-gradient  parameters  defined  by 

.  /  3Uo  3Uo  \ 

\  --  *  (ue  sr  +  sr)  *  LW  "  ueZ>-  (2',6a) 

/  3W  3W  \ 

VN(Ue  JT  *  we  sr)  -  LW  ~  ueZ>  <2-16b> 

and 

N  “  ’  L  =  b  (h^+  K2)  •  (2'17) 


This  new  set  of  equations  (2.10)  to  (2.17)  is  free  of  any  singularities  at 
the  nose  and  moreover  takes  a  well-behaved  form  as  t  •+  0;  thus  they  are 
formally  well -suited  to  Initiate  the  computation  on  thin  prolate  spheroids, 
or  Indeed  any  other  shapes  provided  they  have  paraboloidal  noses.  The 
corresponding  boundary  conditions  are 

u  =  v*w*U  =  V  =  W  =  0  when  y  =  Y  =  0  (2.18a) 

and 

u  -*■  ue,  w  ■*  we,  U  Ue ,  W  -►  Ww  as  y,  Y  «>  (2.18b) 

where  ue,  Wg,  Ue>  Wfi  are  related  by  formulas  equivalent  to  (2.10).  Initial 
conditions  are  also  Imposed  at  the  stagnation  point  SQ,  defined  from  (2.7) 
by 

e  *  0,  e  -  -V0(t)  CV02(t)  +  t2Vg2(t)  tan2a]1/2  i  50*  (2.19) 
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and  are 


U  =  W  =  0 


Y  >  0. 


(2.20) 


In  the  neighborhood  of  this  point, which  we  define  to  be  (XQ ,0)  in  the 
(X,Z)  coordinate  system  we  may  write 

U  =  (X  -  XQ)  ft(Y),  W  =  zft(Y),  V  =  ft(Y)  (2.21) 

and  then  il,  ft  satisfy 

N^2  +  v  a| =  +  v  n5i2  +  v  =  e  2*  +  v  ~t  (2.22) 

N(ft  ♦  ft)  ♦  §  -  0  , 

where 

*r  ••(&)*  (2-23) 

are  evaluated  at  (XQ,0).  These  stagnation  flow  equations  are  a  special 
case  of  Howarth's  equations  (1951). 

While  it  is  true  that  one  may  now  proceed  to  solve  this  new  form  of 
the  three-dimensional  boundary  layer  equations  from  the  forward  stagnation 
point  past  the  nose  and  as  far  downstream  as  there  are  still  some  awkward 
features  about  Initiating  the  calculations;  these  may  be  avoided  by  means  of 
yet  another  transformation.  Since  the  final  form  of  the  equations  Is  so 
convenient  for  the  computer  it  was  decided  to  make  use  of  them  In  spite  of 
the  slightly  more  complicated  algebraic  structure.  The  strategy  underlying 
our  thinking  Is  as  follows.  The  standard- box  method  depends  for  Its  success 
on  knowing  the  solution  on  three  of  the  four  normals  in  which  the  boxes  are 
stacked  at  each  stage.  In  the  formulation  (2.10)-  (2.17)  we  can  determine  the 
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solution  on  Z  =  0  (Cebeci  et  al  1980,  Wang  1970), but  how  can  we  find  the 
solution  on  the  first  normal  of  the  first  station,  say  X  =  0,  Z  =  k,  off 
the  windward  line  of  symmetry?  Of  course.once  we  have  this, the  standard  box 
Is  well  suited  to  finding  the  solution  at  all  the  normals  standing  on  the 
line  Z  =  k,  but  a  special  method  Is  needed  to  get  the  first  solution. 
However.since  we  have  the  solution  on  the  line  of  symmetry  we  could  use  this  pro 
vided  we  changed  to  polar  coordinates  (R*,**)  centered  at  SQ  and 
advanced  the  solution  from  R*  =  0  using  the  box  method  to  carry  us  from 
**  =  0  to  **  =  ir  at  each  new  station  of  R*.  There  Is  another  point. 

Once  we  get  away  from  the  neighborhood  of  the  nose  the  disadvantages  of  the  orig 
inal  system  (2.1)  -  (2.6)  disappear  and  Its  simplicity  both  in  formulation  and 
geometric  Interpretation  are  compelling  reasons  for  basing  our  solution  on  It. 
Thus  we  need  a  transformation  which  converts  the  polar  coordinates  centered 
at  SQ  to  the  'polar  coordinates'  (s+l,e)  centered  at  the  nose  within  an 
acceptably  small  distance  from  It. 

There  are  many  ways  of  achieving  these  aims  and  we  choose  here  to  use 
the  theory  of  coaxal  circles.  We  write 


X  *  x«  +  (n2  "  1)  x„  R(R  +  cos*)a  ,  Z  =  (n2  -  1)  X  Ra  sin*  (2.24a) 


where 


A  =  (1  +  2R  cos*  +  R2)"1 


(2.24b) 


and  n(>l)  is  an  arbitrary  parameter;  later  on  we  shall  set  n  ■  2. 

Other  values  of  n  may  well  be  more  convenient  In  special  circumstances. 
Then  when  R  «  1 


x  -  X  +  XJn2  -  1)  R  cos*  +  0(R2) ,  Z  *  X  In2  -  1)  R  sin*  +  0(R2) ,  (2.25) 


w  =  Ra(1-R2)  si n<t>  Q  +  Ra[cos<|>(1+R2)  +  2R]  T  sin<j> »  (2.27b) 


so  that  when  R  =  1/2.  Q  and  T  sin*  may  be  identified  with  u  and  w. 
Then  the  equations  satisfied  by  Q,  T,  V  are 


(c)  ^-Momentum 


[RQ  Iff  +  s1n*  T  U  +  T2  +  2a(1+R  cos*)QT  -  2aRQ2J 

2 

-  L(S)  Q(WX  -  UZ)  +  V  =  vi{  +  Bj  (R»4> ) . 


Here 


W  si n4>  =  W  and  1-1  sin* 


(2.30) 

(2.31) 


and  g£,  g*  are  pressure  gradient  terms, 'not  reproduced  here,  but  which  may 
easily  be  written  down  in  terms  of  the  main  stream  values  Qe,  Tg  of  Q,  T 
and  these  follow  at  once  from  the  values  of  Ug,  Wg  using  (2127).  The  point 
of  the  relationship  in  (2.31)  is  that  now  the  lines-of-symmetry  equations  can 
be  written  down  at  once  from  (2-28)  -  (2.30)  on  setting  *  *  O.tt  and  remembering 

A  A 

that  W,  1  are  finite.  Thus  the  values  of  Q,T  on  f  *  O.tt  can  be  computed 
separately,  and  in  essence  have  been  (Cebecl  et  al ,  1980),  Independently  of 
the  solution  in  0  <  $  <  it.  Further,  Q  and  T  are  known  at  R  -  0  for 
0  <  <|>  <  ir  from  the  solution  of  the  stagnation  equations  (2.21)  -  (2.23): 


Q(0,*,Y)  -  3Xq  [fif(.Y)  sin2*  -  &(Y)  cos2*] 
T(0,*,Y)  =  3Xq  cos*  [fif(Y)  +  fif(Y)] 


(2.32) 


When  we  couple  known  solutions  at  R  =  0,  and  0  <  *  <  ir,  *  ■  0  and  R  >  0, 

*  *  if  and  R  >  0,  with  the  boundary  conditions  at  y  *  0  and  as  y  •+  «• 
it  Is  clear  that  In  spite  of  the  more  complicated  form  of  (2.28)  -  (2.30)  a 
numerically  attractive  formulation  of  the  original  problem  has  been  achieved. 
We  may  now  use  the  box  method  to  compute  the  solution  systematically 
along  lines  of  constant  R  starting  from  R  =  0  and  increasing  *  from  0  to 
ir  at  each  new  station  of  R.  Further  the  new  form  of  the  equations  is  well- 
suited  for  the  numerical  schemes  proposed  by  Wang  (1974,  1975). 
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Before  embarking  on  the  numerical  solution  It  Is  convenient  to  remove 
u  ,  v  from  the  equations  by  making  appropriate  scale 

oo 

transformations4but  to  save  further  complications  In  notation 
we  shall  simply  set  u,  »  v  *  1  from  now  on. 
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III.  THE  NUMERICAL  PROCEDURE 


The  standard  box  method  for  solving  the  governing  equations  may  be 
Illustrated  by  considering  them  in  their  original  form,  which  we  use  for 
calculating  the  flow  field  when  S  >  2 XQ  or,  more  strictly,  its  equivalent 
in  terms  of  £,  conveniently  defined  as  c  >  ?• 

Difficulties  can  arise  with  this  formulation  when  a  is  small 
because  the  convenient  step  length  interns  of  n  is  much  larger  than  that 
in  terms  of  Y  used  when  5  <  C-  To  avoid  this  problem  Y  may  be  used  as 
variable  until  say  ?  =  0  when  a  further  switch  to  n  could  be  made. 

The  equations  are  reorganized  by  writing 


n  =  y/s 


1/2 


,  s-  j  h, 

A 


and  with  primes  denoting  differentiation  with  respect  to  n» 


u(c,y,e)  =  f(c,n,e),  wU.y.e)  =  gU.n.e) 


(3.1) 


(3.2a) 


f'  *  pU.ri.e)  ,  9'  =  qU.n.e) 


(3.2b) 


*■  _  I_  il  +  l_i£  +  fr1  1 (  \ 
e-h7^7+F~afl  +  f(2i7-ly- 


hi  H  hg  86 


*1 


The  equation  of  continuity  may  now  be  Integrated  to  give 

v  =  /s  (7^-nf-e) 
and  the  momentum  equations  reduce  to 

P'*56,  +Spe-SK2  92  •  s(^  |f) 


(3.2c) 


(3. 2d) 
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(3.2e) 


(3.2f ) 


q'  +  S82  +  sqe  +  sK2  fg 


where  81 ,  e2  are  dimensionless  pressure  gradient  parameters  with  the 
property  that  (3.2e),  (3.2f)  are  automatically  satisfied  in  the  limit 
n  -*•  •  when 


f  -*■  ufi,  g  wg,  p  -*■  0,  q  -*■  0. 
The  additional  boundary  conditions  are  that 


(3.3a) 


f  =  g  =  e  =  0  at  n  =  0 


(3.3b) 


A  three-dimensional  grid  is  now  set  up  in  the  (e,e,n)  space  consisting 
of  straight  lines  in  the  e,e,n  directions  dividing  up  the  domain  of 
integration  into  a  set  of  boxes.  Let  us  conveniently  label  a  particular 
box  by  (i,j,k)  where  the  box  (1,1,1)  has  three  sides  on  the  initial  planes 
S  -  5,  n  =  0  and  e  *  0  respectively.  Further  let  the  e  coordinates  of 
thecornersof  the  (1,j,k)  box  be  where  ^  -  ^.i  =  r^ ,  the  n 

coordinates  be  n^.n^i  where  —  nj-l  =  and  the  6  coordinates 
0n,0n_i  where  0n  -  =  kp  (see  Fig.  2  for  a  typical  box).  It  is 

noted  that  there  is  no  need  for  any  equality  between  any  of  the 


Figure  2.  Finite  difference  molecule  for  the  standard  box. 
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different  r.j,hj,kn.  We  now  suppose  that  f,g,p,q  are  known  at  the  mesh 
points  (c^_i  ,nJ.,en_1),  (ci_1.Hj»en),  (ci  »nj  »en_i )  for  all  j  and  wish  to 
compute  them  at  the  mesh  points  (€^  *nj  »®n) *  The  difference  approximations 
to  (3.2b,c,e,f)  are  written  in  a  standard  form  (see  for  example  Cebeci  & 
Bradshaw  1977)  which  may  be  briefly  described  as  follows.  Equations  (3.2b)  are 
approximated  using  centered  difference-quotients  and  averages  about  the 
mid  point  (?.j  ,nj-yen^ ’  nj-}s  “  (nj-i  +  nj)  •  The  deference  equations  which 
are  to  approximate  (3.2c, e,f)  are  written  about  the  mid-point 
(*J-ynj-y0k-Jp  of  the  cube*  In  centerin9  the  parameters  f,g,p,q  at 
this  point  we  use  formulas  of  the  type 


Vb  *  1  <To ♦  V*  7j  - }  «',*  ♦  *]-'*  *  t"'1  * 


(3.4) 


thus  ascribing  equal  weight  to  the  values  of  f  at  all  corners  of  the  box 

The  centering  of  the  fifth  parameter  e  proceeds  differently.  We  do  not 

tabulate  it  at  the  oomers  of  the  box  but  instead  at  the  mid-points  of  the 

sides  n  =  constant.  Thus  by  analogy  with  (3.4)  we  treat  i.  as 

J 

unknown  and  define 


®j-js  =  ?  (ij +  *j-i}  •  (3-5) 

If  at  any  stage  of  the  computation  we  require  e  on  the  mesh  points,  we  must 
interpolate  among  the  tabulated  values.  This  special  treatment  for 
e  is  necessary  to  avoid  unacceptable  oscillations.  The  derivatives  in 
(3.2c,e,f)  are  approximated  by  formulas  such  as 
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f'*  H7'VVl> 


3f  _  1  Li,n  f1-l,n  .  f1,n-l  «1-l,n-l  ,  -i,n  -i-l,n 

H  '  4?T  [rj  rj  +  Tj  "  fj  +  fj-l  “  fj-l 


+  f 


i  ,n-l 
j-1 


fi-l  ,n-l 
J-1 


(3.6) 


The  difference  equations  together  with  the  boundary  condition  (3.3)  form  a 
non-linear  algebraic  system  for  the  unknown  quantities  ft,n,  g!-’n»  p!’n,  q{,n, 
i.  (all  j)  which  we  linearize  by  Newton's  method  and  solve  the  resulting  linear 

J 

system  using  the  block-elimination  method  discussed  in  Cebeci  and  Bradshaw 
(1977).  In  order  to  start  the  computation* the  values  of  f,g,p,q  but  not  e 
are  required  on  the  planes  9  =  0  and  £  =  £.  These  are  provided  by  the 
windward  line  of  symmetry  solution  already  obtained  as  a  separate  and  self- 
contained  computation  and  by  the  computation  for  the  nose  region.  In  turn 
the  nose  region  can  also  be  solved  by  the  standard  box  provided  we  use  the 
third  formulation  of  the  equations  (2.30)-(3.2)and  is  Initiated  by  the 
stagnation -point  solution  (3.3).  Thus  we  have  a  viable  numerical  procedure 
for  the  Integration  of  the  boundary  layer  equations  which  we  may  use  until 
It  breaks  down  due  to  loss  of  convergence  of  the  iteration  sequence. 

The  seeds  for  this  possibility  are  sown  in  the  arbitrary  assumption  that 
the  integration  must  proceed  in  the  directions  of  5  increasing  and  of  6 
Increasing.  We  note  in  parenthesis  that  the  integration  could  equally  well 
proceed  in  the  direction  of  e  decreasing,  starting  from  the  leeward  line 
of  Integration  since  the  solution  there  can  and  has  also  been  computed 
separately  (Cebeci  et  al  1980),  but  this  does  not  overcome  the  fundamental 
weakness  of  the  standard  box.  For  the  evolution  of  the  boundary  layer  proceeds, 
as  explained  earlier,  by  diffusion  in  the  n  direction  and  by  convection  along 
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the  local  streamline  direction,  at  the  value  of  n  being  considered.  So 
long  as  this  direction  is  roughly  the  same  as  that  of  the  mainstream  and  in  par¬ 
ticular  both  K  and  e  are  Increasing  along  this  streamline,  the  standard  box 
should  be  an  acceptable  procedure;  the  main  feature  that  must  be  checked  is 
that  all  the  local  streamlines  from  points  on  the  lines  e  =  ,  9  *  en 

cross  the  plane  s  =  ^ at  points  inside  the  boxes.  If  not,  the  zone  of 
influence  on  this  line  extends  outside  the  stacked  boxes  but  this  physical 
solution  is  easily  removed  by  decreasing  the  values  of  riskn. 

Suppose  however  the  local  streamline  at  n  =  nj 
is  in  the  direction  of  e  decreasing  while  that  of  the  external  flow 
is  in  the  direction  of  e  increasing.  No  adjustment  of  mesh  lengths  can 
remove  the  unphysical  features  of  the  standard  box  and  a  radical  change  is 
necessary.  We  have  investigated  two  possibilities.  The  first  is  a 
development  of  the  zig-zag  difference  scheme  first  used  by  Cebeci  (1979)  for 
unsteady  flows  and  which  we  shall  refer  to  as  the  zig-zag  box. 

In  this  scheme  we  retain  the  same  procedure  for  solving  (3.2b) .  Suppose 
for  definiteness  we  have  started  the  computation  on  the  windward  line  of 
symmetry.  Then  for  the  other  equations  we  modify  the  algebraic  difference 
equations  only  if  g] <  0.  We  write  (3.2c,e,f)as  algebraic  equations 
centered  at  P,  using  quantities  centered  at  P,Q,R  where 


p=  ^i-ynj-y0n^  »  Q=  ^i’^j-yV^’  R=  ^i-i^j-yV^  <3-7) 


On  the  other  hand  if  the  integration  starts  from  e  =  *  then  the  standard 
box  is  abandoned  locally  only  if  gl’P  >  0.  The  first  situation  occurs  near 
the  surface  n  =  0  and  the  second  in  the  outer  part  of  the  boundary  layer. 

The  pattern  of  the  zig-zag  scheme  in  each  case  may  be  seen  from  Fig  3 
which  shows  a  section  of  the  grid  by  the  plane  n  =  . 

J 
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(a)  (b) 

Figure  3.  Finite-difference  molecule  for  the  zig-zag  box  when  the  solution 
procedure  begins  from  (a)  the  windward  plane,  and  (b)  the  leeward 
plane.  In  each  case  x  denotes  known  values  and  o  computed  by 
the  zig-zag  box. 


To  illustrate  the  centering  of  the  equations  with  this  zig-zag  box 
let  us  consider  a  model  equation  for  (3.2e) 


p'  +  spe  = 


af  g  af 

ar  mF  ao 


(3.8) 


where  we  have  temporarily  added  stars  to  the  metric  functions  to  prevent 
confusion  with  the  step-length  in  the  n  direction.  The  finite  difference 
approximations  to  this  equation  for  the  case  when  the  marching  in  the  e 
direction  is  from  windward  to  leeward  are  (Fig  (3b)): 


p'(P)  +  (spe)(P)  =  (P)  f(P)  (P) 

*  [*„  9<Q>  I?  W)  ♦  I’  -  *„>  9(l»  |f  («)]  <3-9> 
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Here  the  capital  letter  following  each  group  of  numerical  quantities  Indicates 
the  point  at  which  they  are  to  be  centered  and  for  this  purpose  we  use  the 
neighboring  mesh  points  of  e  for  P  and  neighboring  mesh  points  of  e 
for  Q  and  R.  Additionally 


y  -  9n+l  9n 

n  “  Vi  '  V-i 


(3.10) 


Again  we  have  an  algebraic  system  of  equations  for  finding  the  values  of 
the  dependent  variables  at  (£j,nj.en)  (all  j)  which  is  solved  by  the 
procedure  followed  in  the  standard  box  scheme.  A  variation  of  this  method 
has  been  successfully  applied  to  the  computation  of  unsteady  laminar  boundary 
layers  (Cebeci  1979)  and  we  shall  make  use  of  it  here  when  the  cross-flow 
velocity  g  changes  sign  across  the  layer. 


Although  the  zig-zag  box  is  a  definite  improvement  on  the  standard 
box  it  still  contains  some  physical  weaknesses.  Thus  for  the  determination 
of  p  from  (40)  at  (ci  ,0n)  the  values  of  f,  g,  p,  q  at  (^i  ) 

are  not  directly  relevant  since  distributions  are  being  carried  forward  along 
the  local  streamline.  Of  crucial  importance  in  fact  are  the  values  of 
ft  9t  Pt  q,  at  U.j_i  ,0j ,e) ,  this  being  the  point  in  the  plane  s  *  ^  -j 
where  it  intersects  the  projection  of  the  streamline  through  (c^  .Qn ) 
on  the  plane  n  =  n.-. 

J 

Following  the  same  line  of  argument  the  attaching  of  equal  weights  to 
the  centered  mesh  points  on  either  side  of  R  does  not  allow  for  the  possibil¬ 
ity  that  if  this  streamline  passes  close  to  one  of  them  the  other  may  almost  be 
discounted.  For  this  reason  we  use  a  third  form  of  the  Keller-box  method  which 
incorporates  these  physical  arguments  and  which  we  shall  refer  to  as  the 
characteristic  box  (see  Fig.  4).  The  method,  developed  by  Cebeci  and  Stewartson 
(1977),  is  usually  used  for  all  j  when  g  *  0  for  any  j  in  the  calculations 
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starting  from  the  windward  side  (and  vice  versa  from  the  leeward  side)  but  it 
can  be  applied  even  if  g  >0  for  all  j.  The  standard  procedure  is  retained 
for  (3.2b)  and  the  zig-zag  procedure  for  the  continuity  equation  (3.2c).  In 
order  to  obtain  the  difference  approximations  to  (3.2e,f)  we  first  write  them  in 
streamline  form,  that  is: 

P'  +  s81  +  spe  -  sK2g2  =  A  (3.10a) 

q'  +  SB2  +  sqe  +  sKgfg  =  A  |2-  (3.10b) 

where  y  is  measured  along  the  projection  of  the  local  streamline  in  the 
plane  n  =  const  and 

A  =  s  [(f/hf)2  +  (g/h*)2]2  .  (3.11) 

The  value  of  e  is  found  by  drawing  a  straight  line  through  (ci-,nJ-,ek) 

in  a  direction  equal  to  the  mean  of  the  streamline  directions  at  (^_-|» 

n . ,©)  and  (c.,  ,n  •  ,eJ.  This  necessitates  an  additional  measure  of  iteration 
J  *  J  * 

in  the  numerical  scheme.  The  values  of  the  relevant  variables  at  (c^.Oj.e) 
are  found  by  quadratic  interpolation  between  the  known  values  at  (5^_i »nj »®n_i ) » 
(Ci.1,nj,en)  and  .nj.en+1)-  The  equations  (43)  are  then  differenced  at  the 

mid  point  $  between  (c^.nj.e^)  and  (c^_i»nj»0)  (see  Fig. 4),  and  finally  the  value  of 
e  at  »  is  found  by  linear  interpolation  from  its  value  at  P,  where 
equations  to  find  it  already  have  been  written  down,  and  the  values  of  e  at 
smaller  values  of  k.  It  is  believed  that  the  errors  Involved  In  this  procedure 
are  second  order,  i.e.  the  same  as  those  in  the  standard  box  method.  The  method 
has  already  been  used  effectively  In  a  study  of  laminar  and  turbulent  boundary 
layers  on  ship  hulls  (Cebeci,  Chang,  and  Kaups  1978)  as  well  as  three-dimensional 
flows  on  wings. 
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Figure  4.  Finite-difference  pattern  for  the  characteristic  box.  The  heavy 
lines  show  the  pattern  for  the  zig-zag  scheme  for  continuity 
equation,  x  denotes  known  quantities  and  o  those  to  be  computed 
by  the  characteristic  box. 

The  characteristic  box  may  also  be  used  to  extend  the  solution  in  the 
nose  region  away  from  the  line  of  symmetry  in  the  second  form  of  the  equations, 
(2.13)  -  (2. 15). He^a  the  solution  is  known  on  Z  =  0  and  it  enables  us  to  deduce 
the  solution  at  y,  the  first  X-station  of  the  first  grid  lines  adjacent  to  those 
at  Z  =  0.  The  computation  of  e  at  f  is  effected  by  extrapolation  from 
P  using  its  symmetry  properties  on  Z  =  0.  Having  computed  the  solution  at 
this  station  of  X  those  at  other  stations  of  X  follow  using  the  standard 
box.  The  process  is  repeated  at  each  of  the  Z  grid  lines. 

The  characteristic  box  is  complicated  and  requires  a  little  more  programming 
effort  than  either  the  standard  or  the  zig-zag  boxes  but  it  has  a  greater  range 
of  applicability  than  they  or,  we  believe,  any  other  method  available  at 
present.  However  It  does  have  limits;  thus  it  is  likely  to  fall  if  the 
direction  of  v  completely  reverses  in  (43).  This  phenomenon  seems  to  occur 
just  near  the  leeward  portion  of  the  separation  line  i$  and  prevents  us  from 
fixing  its  position  and  properties  as  well  as  we  can  for  the  windward  portion 
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IV.  GENERAL  PROPERTIES  OF  THE  SOLUTION 

The  computations  were  carried  out  for  a  *  3°,  6°,  15°,  30°  and  t  =  1/4. 
For  a  s  6°  the  second  form  of  the  governing  equations  (2.13)  -  (2.15)  was 
used  with  a  switch  to  the  first  form  at  S  *  0.523  (5  =  0.975), but  for  all 
other  values  of  a  the  third  form  was  used  to  Initiate  the  computations. 

The  standard  box  was  used  for  the  Integration  In  the  U.n.e)  variables  until 
the  crossflow  velocity  changed  sign  and  generally  the  characteristic  box  was 
used  subsequently.  A  number  of  studies  were  made  with  the  zig-zag  box  at 
a  =  6°;  it  was  concluded  that  this  variation  is  less  sensitive  to  the 
structural  properties  of  the  boundary  layer  and  discarded.  For  example,  in 
the  integrations  from  the  leeward  line  of  symmetry  a,  the  method  broke  down 
at  £  =  0.44  as  against  £  =  0.35  with  the  characteristic  box.  On  the 
other  hand,  integrations  from  the  windward  line  of  symmetry  *w  broke  down 
at  c  *  0.31  using  either  method.  Various  experiments  were  made  with  the  step 
lengths:  in  the  ^-direction  r.  varied  from  0.01  to  0.1  (Fig.  2),  in  the 
e-direction  kn  varied  from  ZH°  to  10°  (i.e.  0.044  to  0.175)  and  in  the 
n-direction  h.  varied  from  a  minimum  of  0.036  near  the  surface  to  a  maximum 

J 

of  3.7  at  the  outer  edge  when  negative  crossflow  was  well  developed  and  the 

boundary-layer  thickness  had  reached  30.  No  attempt  was  made  to  Investigate 
2 

the  effects  of  h  -extrapolation  or  other  deferred  approaches  to  the  limit  on 
the  results  obtained. 

The  velocity  profiles  Individually  show  no  new  general  features  from  those 
reported  by  earlier  authors  (Wang  1974a, b,c,  1975,  Gelssler  1975)  and  will  not 
be  discussed  in  detail  here.  The  growth  of  the  boundary  layer  In  the  region 
of  crossflow  reversal  has  already  been  noted,  but  in  contrast  to  the  results 
obtained  by  Cebecl  et  al  (1980)  for  the  paraboloid  (t  =  0)  we  find  that  the 
thickness  of  that  part  of  the  boundary  layer  In  which  g  <  0  does  not 
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Increase  significantly.  Thus  when  a  *  30°  that  part  of  the  boundary  layer 
never  extended  beyond  n  *  2  even  though  uniform  conditions  might  not  be 
achieved  until  n  *  30. 

The  discussion  is  concentrated  on  the  principal  features  of  the  boundary 
layer,  namely  the  skin  friction  and  the  effective  blowing  velocity  which  it 
induces  on  the  external  flow.  From  the  equation  of  continuity  (2.1)  we  know 
that  v  asymptotes  to  a  linear  function  of  y  as  y  -*■  «  on  the  boundary 
layer  scale.  Let  us  define 


vje.e)  =  Lt 

y  -►  oo 


v  + 


h^ 


3x  ^h2U^ 


3 

ase 


(h-|W) 


(4.1) 


Then  the  perturbation  to  the  inviscld  flow  caused  by  the  boundary  layer  may 
be  regarded  as  a  replacement  of  the  condition  of  zero  normal  velocity  on  the 
body  by  the  requirement  that  it  be  equal  to  v^.  Equivalently,  the  condition 
that  the  body  be  a  streamline  surface  is  changed  to  the  condition  that  the 
streamline  surface  through  the  forward  stagnation  point  SQ  is  at  a  distance 
6*  from  the  body  where 


I*  <h2V*>  +  ale  <hlV*>  =  hlh2v-  i  <4-2> 

6*  can  be  found  by  solving  this  equation  with  the  initial  condition 

6*  =  0  at  SQ.  By  reference  to  (3.2)  we  may  connect  v  with  e  and  hence, 

after  setting  uw  =  a  =  v  =  1  and 


we  have 

v.  =  E.s*  .  (4.4) 


The  principal  quantities  on  which  we  shall  base  our  discussion  of  the  solution 
properties  are  pw,  qw  and  E^,  the  components  of  the  skin  friction  in  the 
5,e  directions  being 


(4.5) 


s-i*lpU,n,0),  q(C,n,0)]  =  s’*(Pw»<Vf)  * 

The  equivalent  physical  quantities  may  now  easily  be  deduced  on  restoring 
the  appropriate  values  of  um,  a,  v. 

In  presenting  the  results,  we  shall  pay  closest  attention  to  the  solu¬ 
tions  for  o  =  6°  and  for  a  =  30°,  the  first  because  the  ok  is  well -developed 
yet  a  laraelv  comolete  solution  can  be  found, and  the  second  because  the  ok  has 
penetrated  almost  to  the  nose  (5  “  -0.83)  and  furthermore  the  accessibility 
boundary  differs  most  clearly  from  on  the  leeside.  The  next  sig¬ 
nificant  stage  in  the  development  of  1 ^  with  increasing  o  occurs  at  o  <*  42°, 
where  joins  up  with  an  incipient  nose  separation  and  the  whole  of  the 
leeside  of  becomes  inaccessible.  Careful  studies  of  the  boundary  layer 
have  been  made  in  this  situation  by  Wang  (1974b)  and  we  see  no  reason  to 
repeat  them. 

In  Figs.  5,  6  we  show  the  variation  of  pw,  qw,  E^  over  the  majority  of 
the  accessible  region  lying  upstream  of  «s.  Perhaps  the  most  interesting 
result  is  that  <  0  over  a  significant  part  of  the  region  of  negative 
crossflow.  Thus  the  commonly  held  view  that  boundary  layers  under  adverse  pressure 
gradients  act  as  Injection  sources  with  respect  to  the  external  flow  is  seen 
not  to  be  universally  correct  and  fails,  in  particular,  near  the  leeside  line 
of  symmetry  a  when  the  pressure  gradient  is  as  we  have  prescribed.  It 
does  not  necessarily  follow  of  course  that  6*  <  0  in  this  region, but  we  note 
that  for  the  paraboloid  (t  =  0)  discussed  by  Cebecl  et  al  (1980)  the  value 
of  <5*  tends  exponentially  to  -»  at  Infinite  distances  downstream  from  the 
nose  according  to  the  asymptotic  formula  derived  there..  The  negative  values 
of  E^  may  be  Inferred  from  Wang's  studies  (1974a ,p. 38)  when  it  Is  borne 
in  mind  that  in  his  notation 
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Figure  5a.  Variation  of  the  longitudinal  component  of  wall  shear  parameter  pw  around  the  circumference 
of  the  body  for  different  c-values,  a  =  6°. 
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Figure  5b.  Variation  of  the  transverse  component  of  wall  shear  qw  around  the  clrcumfe 
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Figure  6b.  Variation  of  the  transverse  component  o'  wall  shear 
qw  around  the  circumference  of  the  body  for  differei 
5-values,  a  =  30°. 


(4.6) 
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The  principal  Interest  In  the  plots  of  the  skin-friction  lines  is  In 
the  location  of  i  and  its  nature.  The  location  of  the  line  of  zero 
crossflow  skin  friction  is  shown  in  Fig.  7  for  the  various  values  of  a. 

Good  agreement  is  obtained  with  the  specific  results  of  Geissler  (1975)  at 
ot  =  15°  and  of  Patel  and  Choi  (1979)  at  a  =  6°,  but  there  are  significant 
differences  with  Wang's  data  (e.g.  1975). 

In  Fig.  7  we  also  plot  our  best  estimates  of  the  separation  lines 
for  the  same  values  of  a.  The  principal  properties  of  the  solution  near  are 

discussed  in  detail  in  the  next  two  sections  but  broadly  on  the  windward  side  its 
shape  is  deduced  by  a  combination  of  numerical  analytic  arguments  but  on  the 
leeside  it  is  essentially  determined  by  the  failure  of  the  numerical  scheme 
as  e  decreases.  In  addition  we  have  superposed  on  these  diagrams  the 
external  streamlines  and  we  see  that  at  ot  =  15°,  30°  and  on  the  leeside  of 
they  are  pointing  into  the  region  of  computation.  In  accordance  with 
Raetz's  notion  of  zones  of  influence,  it  follows  that  the  solution  in  part  of 
this  region  cannot  be  completed  with  the  information  supplied  in  this  report. 

For  information  is  also  necessary  from  the  separated  region  downstream  of 
which  we  are  unable  to  compute.  Thus  the  definition  of  separation  as  the 
limit  of  the  accessible  region  from  the  nose  is  apparently  not  compatible  either 
with  it  being  a  limiting  streamline  or  an  envelope  of  streamlines.  We  shall 
return  to  this  point  below. 
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Figure  7.  (b)  a  =  6°. 
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4.1  The  Separation  Line  When  a  - 


An  important  goal  of  our  numerical  studies  Is  to  clarify  further  the  nature 

of  the  separation  line  a$.  For  this  reason  particular  attention  was  paid  to  the 

properties  of  the  solution  in  Its  neighborhood  when  a  =  6°.  The  conclusions 

we  are  able  to  come  to  for  this  special  angle  of  incidence  may  then  be  used 

to  comment  on  the  nature  of  as  for  other  values  of  a.  There  is  no  doubt 

that  when  a_  intersects  either  the  windward  line  of  symmetry  a  (at 
s  w 

S  =  0.84)  or  the  leeward  line  of  symmetry  a  (at  c  “  0.72),  the  solution 
develops  a  singularity  of  the  Goldstein  type  (Goldstein  1948,  Brown  1965). 

We  may  also  expect  (Buckmaster  1971)  that  at  these  points  a$  is  perpendicular 
to  aw,a  and  further  that  in  their  neighborhoods  it  coincides  with  the 
accessibility  boundary  a^. 

There  are  several  conflicting  views  about  the  flow  properties  at  general 
points  of  a  .  Let  us  suppose  first  of  all  that  A  is  a  continuous  curve 
joining  U,e)  *  (0.84,0)  to  (0.72, it)  but  is  not  necessarily  monotonic.  Thus 
we  postpone  discussion  of  Wang's  (1976a)  concept  of  open  separation  until 
later.  In  this  case  there  has  been  general  support  for  Hayes  (1951)  concept  of  a$, 
or  strictly,  the  normals  to  the  surface  through  i$,  as  a  bound  of  those 
points  for  which  the  boundary  layer  can  be  computed  starting  from  the  forward 
stagnation  point  and  making  use  only  of  conditions  at  n  =  0  and  as  n  -*■  «. 

Beyond  this  bound  the  solution  depends  in  some  way  on  conditions  further 
downstream  and  In  particularly  at  the  rear  stagnation  point.  We  concur  with 
the  notion  that  our  forward  Integration  procedure  cannot  continue  the  solution 
beyond  a$  but  we  have  strong  reservations  about  the  implication  that  It  may 
be  continued  as  far  as  a$  In  all  cases.  When  a  £6°,  however, the  Identification 
of  a A  and  a$  Is  In  our  view  largely  complete. 

The  main  controversy  about  Ag  In  the  literature  centers  on  whether  It 
is  an  envelope  of  limiting  streamlines, i.e.  skin-friction  lines,  (Eichelbrenner 
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and  Oudart  1954,  1973,  Maskell  1955)  or  is  itself  a  limiting  streamline 
(Lighthlll  1963).  The  distinction  is  not  academic.  In  fact  It  is  of  crucial 
Importance  in  the  study  of  practical  boundary  layers,  whether  laminar  or  tur¬ 
bulent,  because  if  the  first  view  is  correct  then,  unless  special  preventa¬ 
tive  measures  are  taken,  a  singularity  must  develop  in  the  solution  at  all 
points  of  In  this  event  the  slope  of  the  displacement  surface  also 
becomes  Infinite  at  and  the  boundary  layer  exerts  a  serious  modifying 
influence  on  the  mainstream.  Thus  the  hierarchical  method  of  studying 
high  Reynolds  number  flows  fails,  just  as  it  does  in  two  dimensions. 
Experimentally  Inspired  statements  about  are  therefore  not  germane  to 
the  question, for  In  a  practical  flow  the  solution  must  be  regular  and  Lighthill' 
notion  of  it  being  a  limiting  streamline  Is  correct.  However,  it  is  at 
present  virtually  Impossible  to  prescribe  the  correct  pressure  gradient  and, 
if  it  turns  out  that  a  general  pressure  gradient  leads  to  a  singularity,  there 
are  Important  consequences.  For  example,  a  general  interaction  theory  of  the 
kind  Investigated  by  Smith  0977,  1978)  will  almost  certainly  be  necessary  to 
treat  the  solution  in  the  neighborhood  of  sl$. 

We  can  see  how  the  singularity  must  develop  at  *s  if  it  is  an  envelope 
of  limiting  streamlines  as  follows.  Suppose  that  is  given  by  ?  = 

for  all  e  and  near  £  =  z  the  two  components  of  skin  friction  are 

t5  -  Us  -  C)6  te  =  1  6  >  0  (4.7) 

Then  the  limiting  streamlines  are  the  curves 

6  "  err^T  (cs  ”  e)1'*  (4.8) 

where  e-|  Is  a  parametric  constant.  Thus  If  8  <  1  these  curves  touch 

i%  at  e  *  0^  »  whereas  If  3  *  1  (corresponding  to  a  regular  solution) 

they  meet  It  at  e  =  ». 

After  careful  examination  of  the  numerical  data  near  the  termination  of 
our  calculations,  we  obtained  strong  but  not  totally  convincing  evidence  that 


the  separation  line  Is  an  envelope  and  we  are  led  to  make  the  hypothesis  that 
this  property  holds  over  the  whole  of  Its  accessible  length.  The  notion  Is, 
however,  not  without  Its  difficulties,  especially  on  the  leeward  side. 

Let  us  begin  by  looking  at  the  calculations  in  the  neighborhood  of 
on  the  windward  side,  i.e.  for  e  <  110°  (Fig. 7b).  As  explained  in  Section  III, 
when  we  integrate  the  equations  at  a  given  value  of  £  for  increasing  values 
of  e  starting  from  tw,  where  e  =  0,  our  procedure  is  to  use  the  standard 
box  until  the  crossflow  velocity  w(=g)  changes  sign.  Then  we  switch  either 
to  the  zig-zag  box  or,  more  usually,  to  the  characteristic  box,  both  of 
which  take  into  account  to  an  increasing  degree  the  fact  that  information 
is  being  fed  to  the  new  (£,e)  station  from  larger  as  well  as  smaller  values 
of  e.  As  e  increases  further,  a  stage  is  reached  when  neither  of  these 
boxes  can  be  used  through  lack  of  information  from  the  next  e  station  at 
the  previous  £  station  (the  points  (i-1,  j,  m+1)  in  Fig.  3b)  or  because 
P  becomes  negative  or  because  the  iterations  fail  to  converge.  In  the 
first  eventuality  the  standard  box  is  used  for  one  further  e  station  and 
the  calculations  at  this  £  station  are  terminated.  This  procedure  is 
rather  crude  and  it  is  desirable  for  a  full  understanding  of  the  separation 
singularity  that  a  superior  way  of  overcoming  the  difficulty  be  found. 

A  new  set  of  values  of  p  ,  q  ,  E  is  displayed  in  Table  1;  they 

W  W  ® 

correspond  to  a  representative  value  0.7  of  £  and  at  values  of  e  near  the 
termination  of  the  calculations  which  occurred  at  e  =  tt/3  (*60®).  At  this 
station  due  to  the  lack  of  any  information  about  the  solution  at  £  =  0.69 
and  e  =  62.5°;  the  wall  shear  parameter  p  became  negative  after  only 

W 

two  iterations  causing  the  iterations  to  diverge.  For  this  data  we  infer  that 
the  separation  line  passes  close  by  (0.7,  n/3)  and  that  the  solution  is  not 
smooth  there. 
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Table  1.  Computed  values  of  p  ,  qw,  E_  at  £  =0.7, a  =  S'*.  The  e  step 
Length  in  radians  is  k  *  0.O437  and  the  values  of  E«  were 
actually  calculated  at  midway  stations  In  both  c  and  e. 

e°  P„  q,  E.  (P»tK>2  I  pw  ~  qw  E-2 

42.5°  0.413  +0.002  2.091  0.172  0.136  0.229 

45°  0.394  -0.005  2.412  0.154  0.136  0.172 

47.5°  0.369  -0.014  2.857  0.136  0.137  0.123 

50°  0.343  -0.024  3.468  0.118  0.139  0.083 

52.5°  0.313  -0.035  4.254  0.091  0.139  0.055 

55°  0.279  -0.049  5.415  0.069  0.140  0.034 

57.5°  0.233  -0.066  7.692  0.045  0.144  0.016 

fin0  0.104  -0.098  16.94  0.005  0.133  0.003 


I  Pw  "  qw 

E'2 

00 

0.136 

0.229 

0.136 

0.172 

0.137 

0.123 

0.139 

0.083 

0.139 
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We  now  test  the  hypothesis  that  the  separation  line  Is  an  envelope  of 
limiting  streamlines.  On  this  assumption,  Brown  (1965)  has  studied  the 
nature  of  the  solution  near  separation  and  her  conclusions  may  be  interpreted 
in  the  present  context  as  follows:  Let  (Cs*es)  be  a  point  P$  on  the 
separation  line  which  makes  an  angle  with  the  ^-direction. 

At  f*s  the  component  t$  of  the  skin-friction  along  the  direction  of  i 
finite  while  the  normal  component  Tn  vanishes.  Her  theory  predicts  that 
at  a  point  Q  on  the  normal  to  A  through  P$ 


while 


xs  a  A$  +  Bs(QPs)“  +  ...  when  QP$  «  1 


=  c  fQPJH  +  ...  when  QP,  «  1 


(4.9a) 

(4.9b) 


where  A.  Be,  C.  are  numbers  with  C  >  0.  We  now  tabulate  p  -  tane 

S  S  S  5  w  w 

and  p  tane  -q  as  functions  of  9  for  various  e.  These  functions  are  not 
w  w 

quite  the  in,  t$  of  (4.9)  since  the  e-direction  is  not  normal  to  but  these 

differences  are  not  thought  to  be  significant  and  the  final  value  of  e  chosen 

must  be  consistent  with  the  predictions  of  the  separation  point  at  5  =  0.70 

and  at  neighboring  values  of  £.  We  conclude  that  the  best  value  of  tane  ®  -1/3, 

2 

corresponding  to  ©s  t  -19°  and  display  in  Table  1  the  values  of  (pw  +  1/3  qw) 
and  of  1/3  pw  -  qw-  It  is  seen  that  (pw  +  1/3  qj2  is  roughly  linear,  in 
agreement  with  (1),  except  very  near  e  =  60°  and  we  infer  that  separation 
occurs  at  0S  =  60.5°.  A  more  cautious  conclusion  is  that  (pw  *  1/3  qj<\-  (e$  -  e) 
where  1/3  ±  n  <  1/2  as  9+  es~,  but  in  view  of  the  complexity  of  the  numer¬ 
ical  program,  the  slight  irregularity  visible  in  the  data  of  Table  1  and  the 
difficulty  in  reproducing  singularities  by  numerical  integration,  we  are  broadly 
satisfied  that  these  results  are  consistent  with  Brown's  theory.  The  function 

-  1/3  qw  is  almost  constant  except  right  near  e  «*  60°,  where  the  second  term 
w  w 

ir.  (4.9a)  might  be  important  but  it  is  more  likely  that  the  program  Is  inadequate 

_2 

to  predict  this  term.  Finally  we  display  E^  .which  should  be  linear  in 
e  with  a  zero  at  e  =  e$.  This  is  reasonably  the  case  except  near  e  =  e$ , 
but  the  predicted  value  of  eg  from  this  set  of  values  in  Table  1  is  approx¬ 
imately  59°, rather  less  than  that  predicted  from  p  ,  q  ,  especially  when  it  is 

W  W 

borne  in  mind  that  E  is  tabulated  for  5  =  0.695. 

09 

After  repeating  these  arguments  at  all  the  5-stations  on  the  windward 
side  we  are  able  to  determine  *s  from  5  =  0.31, where  es  *  113?  to  t  «  0.83, 
where  e$  =  0,  and  confirm  that  its  shape  is  consistent  with  the  choice  of 
e  in  tables  similar  to  Table  1  .  We  conclude  that  on  the  windward  side 
is  is  Indeed  an  envelope  of  limiting  streamlines  and  that  the  skin  friction 
T.  on  It  (=p  sin©  -  q  cos©)  varies  from  »0.13  at  ?  8  0.31  to  a  maximum 
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»0. 15  at  £  =  0.6,  decreasing  afterwards  to  zero  when  *s  Intersects  the 
windward  line  of  symmetry.  The  shape  of  the  separation  line  is  shown  in  Fig. 7b 
and  the  variation  of  along  It  in  Fig.  8.  Along  fcw»  Is 
perpendicular  to  in  agreement  with  Buckmaster's  (1971)  theory  of  the 
structure  of  the  stparation  line  in  this  neighborhood. 

On  the  leeward  side  the  situation  is  much  more  difficult  to  Interpret 
with  confidence.  First  of  all,  in  the  integration  from  e  *  *  no  separation 
is  observed  for  £  <  0.35,  even  though  from  the  windward  side  it  occurs  at 
£  *  0.31  and  when  0.31  <  £  <  0.35  the  two  integrations  give  results  closely 

in  agreement  for  e  <  100°.  If  separation  does  occur,  it  Is  very  weak  and  may  have 

been  missed  in  the  numerical  calculation  from  i  because,  even  with  the  use  of 
the  characteristic  box  there  is  a  slight  loss  of  precision  in  integrating  against 
the  cross-stream  direction  outside  the  boundary  layer,  especially  as  g  >  0 
for  all  0  when  n  >  3  and  5  <^0.35.  The  first  onset  of  separation,  at 

£  =  0.35,  occurs  when  0  <  110°  which  is  less  than  the  value  of  es 

at  £  =  0.31  according  to  the  calculation  from  the  windward  side.  We 
believe  that  the  discrepancy  between  the  windward  and  leeward  integrations 
is  largely  the  fault  of  the  latter  and  that,  bearing  in  mind  the  very  sharp  peak 
in  at  £  *  0.30  (see  Fig.  5c),  separation  occurs  first  close  to 
(0.31,  113°). 

For  5  ^0.35  the  integration  from  0  =  tt  breaks  down  at  values  of 
0  >_  110°  in  a  curious  way.  As  0  decreases  pw  decreases  slowly, 
q  increases  steadily  and  E  remains  almost  constant.  As  the 
calculations  draw  to  an  end  pw  begins  to  fall  more  rapidly  towards 
negative  values,  sometimes  sharply,  and  the  Integration  then  terminates 
either  by  a  failure  to  converge  or  a  lack  of  information  from  the  previous 


C  station,  even  within  the  requirements  of  the  standard  box.  These  causes  of 
failure  occurred  roughly  in  the  ratio  2:3,  respectively.  We  give  a  set  of 
results  in  Table  2  for  5  =  0.6  typical  of  the  situation  where  termination 
occurs  through  lack  of  the  necessary  data.  The  sudden  onset  of  breakdown  makes 
unprofitable  a  study  of  the  solution  properties  near  of  a  similar  kind  to 
that  just  carried  out  for  the  windward  side.  It  is  clear  that  something  odd  is 
happening,  but  Brown's  theory  is  of  little  direct  help. 

Table  2.  Computed  values  of  p  ,  q  ,  E„  at  £  =0.6,  a  =  6°.  The  step 
length  in  radians  is  wkn  -  0.0437  and  the  values  of  E«o  are 
actually  calculated  at  midway  stations  in  both  £  and  e. 


e 

Pw 

«w 

E 

OO 

165 

0.156 

-0.120 

-6.94 

162.5 

0.154 

-0.139 

-6.84 

160 

0.151 

-0.158 

-6.76 

157.5 

0.148 

-0.176 

-6.67 

155 

0.146 

-0.193 

-6.52 

159.5 

0.144 

-0.210 

-6.36 

150 

0.121 

-0.226 

-6.08 

147.5 

0.120 

-0.215 

1 

-4.99 

Let  us  suppose  that  is  here  also  an  envelope  of  limiting  streamlines 
the  alternative,  that  It  Is  Itself  a  limiting  streamline  seems  excluded  by  the 
breakdown  of  the  iterative  process  at  the  final  value  of  e.  Then  either  the 
limiting  streamlines  turn  upwards  before  touching  as  In  Fig.  9a,  or 
turn  back,  as  in  Fig  .  9b).  The  first  alternative  seems  unlikely 
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Figure  9.  Two  optional  limiting  streamline  patterns  near  leeside  part  of 

i  . 
s 

from  the  available  data  since  q..  shows  no  signs  of  becoming  positive.  The 
second  requires  pw  to  change  sign, which  is  more  consistent  with  the  data  as 

the  final  computed  value  of  p  is  often  less  than  half  its  value  near  e  =  it. 

It  also  agrees  with  the  Buckmaster  theory ,  which  applies  when  intersects  &. 

If  the  Brown  structure  holds  very  near  to  the  skin  friction  x$  along 
ts  is  closely  given  by  the  last  reasonable  values  of  pw,  qw  and  the  slope 
of  ts  may  be  inferred  from  the  line  of  breakdown  of  solutions.  We  can  also  estimate 
the  variation  of  r$  over  the  range  0.35  <  c  <  0.71  on  the  leeward  side;  it  is 
displayed  In  Fig.  8.  The  position. of  *s  may  be  inferred  as  well 

within  an  accuracy  of  about  1°.  A  necessary  requirement  for  obtaining  the 

solution  right  up  to  Is  that  the  local  streamlines  on  the  normals  to 
the  surface  through  *s  are  all  directed  out  of  the  origin  of  integration. 

The  most  likely  streamlines  to  violate  this  requirement  are  those  external  to 
the  boundary  layer,.but  for  a  =  6°  our  estimate  for  ns  does  not  lead  to  a 
contradiction  due  to  this  cause  except  near  the  ok.  The  external  streamlines  are 
also  indicated  in  Fig.  7b  and  are  pointing  out  of  the  computed  region  at  the 
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estimated  position  of  ts-  We  infer  that  for  o  =  6°  the  boundary  of  the 
region  accessible  from  the  forward  stagnation  point  is  largely  defined  by  the 
position  of 

For  the  shape  of  near  the  ok  at  £  =  0.31  three  possibilities 
were  considered  and  are  shown  in  Fig.  10.  Our  preference  is  Fig.  10a  and  if  we 
are  right,  the  complicated  structures  of  the  limiting  streamlines  near  the  ok 
makes  their  computation  very  difficult.  Some  comments  are  in  order  about  (b) 
and  (c).  The  option  (b)  is  inspired  by  the  notion  of  "open  separation" 
which  Wang  (1974b)  introduced  in  his  study  of  the  solution  when  a  =  30° 
and,  if  it  occurs.it  implies  that  t  vanishes  at  some  point 


Figure  10.  Three  optional  limiting  streamline  patterns  near  the  Ok. 
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of  i%.  The  option  (c)  Is  a  weak  boundary  layer  collision  (shown  by  the 
dotted  line)  and  would  be  Identified  in  our  calculations  by  a  separation 
line  overlap  when  computed  from  i  and  *w  for  fixed  5.  It  has  been 
found  in  boundary  layer  studies  on  rotating  spheres  (Banks  1 976)  and  also 
in  the  entrance  region  of  a  curved  duct  (Stewartson,  Cebecl,  and  Chang, 

1980)  ,  where  it  is  preceded  by  a  singularity  in  the  displacement 
thickness.  Our  calculations  (Fig. 5c)  indicate  that  a  singularity  in  En  is 
very  close  to  (0.30,  113°),  but  they  are  not  as  smooth  in  this  neighborhood 
as  we  should  like  and  a  clarification  of  the  structure  might  well  follow 
from  a  refinement  in  our  numerical  work. 

4.2  The  Separation  Lines  for  a  =  3°,  15°,  30° 

When  a  =  3°  the  first  signs  of  irregularity  are  the  rapidly  increasing 
values  of  E^  for  e  %  105°  as  £  approaches  0.545.  The  windward  integration 
breaks  down  at  (0.545,  110°)  with  pw  rapidly  approaching  zero  as  this  point 
is  reached.  The  leeslde  computation  may  be  continued  to  0  =  0°  at  e  =  0.545 
but  the  maximum  value  of  is  11  at  0  =  100°  suggesting  that  separation 

is  near.  At  e  =  0.60  the  leeside  computation  breaks  down  at  0  *  120°  and 
at  5=  0.625  no  solution  can  be  obtained  since  separation  has  already 
occurred  at  0s  180°.  The  windward  calculations  behave  similarly  to  those 
for  a  *  6°.  We  Infer  that  ts  is  closed,  and  the  normals  through  it  very  largely, 
if  not  entirely,  define  the  accessibility  limit  1^.  It  is  noted  that  the 
general  shape  of  was  forecast  by  Wang  (1976). 

When  a  =  15°  the  first  sign  of  difficulty  in  the  numerical  computation 
occurs  at  c  =  -0.475  when  an  irregularity  develops  in  the  windward 
integration  at  0  =  135°;  in  addition  Ew(c,135)  reaches  a  value  of  14  at 
this  point, which  suggests  that  separation  is  near.  Breakdown  first  occurs 
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at  (-0.425,  135°)  as  the  calculations  advance  from  the  windward  side.  At 
£  =  -0.400  the  integrations  in  both  directions  break  down,  at  130°  and  135°,  the 
leeside  one  because  the  iterations  exceed  the  maximum  allowed.  For  -0.4  <  £  <  0.1 
the  leeside  computations  come  to  an  end  at  e  =  140°,  and  at  £  =  -0.1  the 
solution  at  e  =  140  becomes  unacceptable  because  |q|  >  0.1  at  the  outer 
edge  of  the  boundary  layer  in  a  dramatic  fashion,  no  such  phenomenon  being 
observed  at  0  =  145°.  Thereafter  the  lower  limit  of  0  slowly  increases  with 

£  reaching  165°  at  £  =  0.2.  The  physical  explanation  is  that  the  integra¬ 

tion  at  the  lower  limit  of  0  should  make  use  of  information  from  even  smaller 
values  of  0  in  the  outer  part  of  the  boundary  layer  and  cannot  do  so  because 

this  information  has  to  come  from  the  separated  part.  From  Fig.  7c  the  true 

lower  limit  of  e  defining  the  accessible  region  is  almost  certainly  given 
by  the  external  streamline  from  the  ok,  and  the  position  of  on  the 
leeward  side  is  irrelevant  to  the  computations.  The  computations  were 
arbitrarily  terminated  at  £  =  0.200;  although  the  boundary  layer  is 
thick  (^18)  there  is  no  new  fundamental  reason  why  they  cannot  be 
continued. 

On  the  windward  side  there  are  no  new  features  of  the  solution  that  need 
special  comment;  it  is  interesting  that  *s  nearly  coincides  with  the  line  of 
zero  qw>  The  ok  region  of  is  difficult  to  be  definite  about,  there  being 
again  the  three  possibilities  of  Fig.  10  to  consider.  In  this  case  also  our 
inclination  is  to  favor  Fig.  10a,  that  is  closed  in  this  neighborhood 
but  the  evidence  is  very  weak.  Wang  (1976)  favors  open  separation  at 
a  =  30°  over  a  substantial  proportion  if  not  all  of  the  leeward  part  of 
In  our  view  the  notion  of  inaccessibility  preempts  a  definite  decision  on 
this  question. 
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When  a  =  30  ,  the  first  sign  of  difficulty  occurs  at  (-0.85,  155°) 
when  reaches  a  peak  of  =19  suggesting  that  separation  is  near.  From 
the  windward  side  breakdown  first  occurs  at  c=  -0.83  and  e  =  140°  where 
Pw  <  0.  On  the  leeside  breakdown  also  occurs  at  5=  -0.83  but  at  3=  146°; 
these  points  are  sufficiently  close  together  that  their  actual  coincidence  may 
reasonably  be  inferred.  Continuing  with  the  leeside  solution,  the  ok  shows  some 

indication  of  being  similar  to  that  for  a  =  6°  and  then  the  breakdown  value 
of  e  slowly  increases  with  e  reaching  156°  at  c  =  -0.74  at  which  point 
| q |  >0.1  at  the  outer  edge  of  the  boundary  layer.  As  for  o  =  15°,  the 
reason  is  that  the  boundary  of  the  accessible  region  from  stagnation  is  prob¬ 
ably  upstream  of  £$.  We  note  that  the  breakdown  in  the  solution  for 
-0.83  <  5  <  -0.74  on  the  leeside  cannot  be  directly  associated  with  is  and 
also  happens  to  occur  at  points  practically  coincident  with  the  external  stream¬ 
line  through  the  ok.  The  shapes  of  this  external  streamline  of  £$  on  the 
windward  side,  and  of  the  breakdown  curve  on  the  leeward  side  are  displayed  in 
Fig.  7d  to  demonstrate  this  property.  On  the  windward  side  the  behavior  of  *s 
is  similar  to  that  for  other  values  of  a  and  we  shall  not  comment  on  it. 


V.  DISCUSSION 


The  generalization  of  the  Keller  box  method  of  solving  parabolic 
differential  equations,  to  take  account  of  the  varying  direction  of  the 
local  streamlines  from  that  of  the  limiting  streamline  at  the  body  to  the 
external  streamline  Including  any  overshoots,  by  the  concept  of  the  characteristic 
box  has  enabled  us  to  continue  accurate  solutions  of  the  boundary  layer 
equations  further  than  has  been  possible  hitherto.  In  fact  we  believe  that 
the  method  can  be  used  virtually  up  to  the  boundary  of  the 
accessible  region  from  0. 

At  zero  incidence  and  are  coincident  circles  defined  by 
£  ^0.79  (Wang  1970).  At  a  =  3°  the  two  are  still  coincident  but  are  no 
longer  circles  (Fig.  7a)  being  moved  forward  on  the  leeside  to  £  ~  0-60 
(i.e.  for  e  >  100°)  and  curving  back  to  £  ;v0.76  at  e  =  0.  A  large 
region  of  weak  reversed  circumferential  flow  develops  extending  almost  to 
£  =  0  near  i.  There  is  little  sign  of  an  ok  developing  at  this 
stage,  but  when  a  reaches  6°  it  is  quite  pronounced, extending 
back  to  £  *  0.31  at  e  =  110°.  There  is  some  uncertainty  at  the  tip 
of  the  ok  but  elsewhere  we  are  confident  that  the  separation  line  defines  the 
limit  of  the  region  of  accessibility  since  the  streamlines  on  the 
normals  to  either  touch  this  surface  or  are  directed  out  of  It 
I.e.  away  from  0.  The  region  of  cross-flow  reversal  extends  further  towards 
the  nose  and  the  phenomenon  becomes  stronger  so  that  the  characteristic  box 
Is  essential  If  the  whole  region  upstream  of  Is  to  be  computed.  On  the  wind¬ 
ward  side  cross-flow  reversal  Is  soon  followed  by  separation  and  a  careful 
study  of  the  natures  of  the  solution  near  supports  the  view  that  it  is 
singular  at  *s  as  described  by  Brown.  This  Irregularity  Is  best  seen  in  the 
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behavior  of  E^  (Fig.  5c)  which  is  clearly  infinite  at  separation.  On  the 
leeward  side  the  structure  is  more  obscure  but  the  limiting  streamlines 
appear  to  have  to  turn  until  they  are  almost  pointing  in  the  opposite 
direction  to  those  of  external  streamlines  above  them.  The  structure  of 
the  neighborhood  of  the  ok  is  very  certain  but  what  evidence  there  is  suggests 
that  is  smooth  and  continuous  there.  At  the  ok  the  direction  of  the 
external  streamline  is  unfavorable  on  the  leeward  side  and  so  according  to 
Fig.  7b  there  may  well  be  a  small  length  of  £^  not  coincident  with  £$, 
which  is  then  uncomputable. 

At  a  *  15°  and  30°  the  windward  part  of  £$  is  very  close  to  the 
cross-flow  reversal  but  there  seems  little  doubt  that  £s  is  an  envelope 
of  limiting  streamlines, the  most  convincing  evidence  coming  from  the 
behavior  of  E^.  The  singularity  In  E^  is  strongest  at  the  tip  of  the 
ok  and  on  the  leeward  side  the  solution  breaks  down  almost  exactly  on  the 
external  streamline  through  this  point.  No  difficulty  at  all  is 
experienced  in  integrating  the  equations  from  £  provided  the  characteristic 
box  is  used  so  long  as  we  are  upstream  of  this  line.  Here  however,  and 
quite  suddenly, convergence  of  p,q  Is  lost  at  the  outer  edge  of  the  boundary 
layer  as  would  be  expected  since  the  corresponding  external  streamline  is 
carrying  the  data  of  Its  numerical  singularity  at  the  ok  with  it.  Thus  the 
limit  £^  of  the  accessible  region  on  the  leeward  side  Is  not  provided  by 
£$  but  by  the  external  streamline  through  the  ok.  Not  only  have  we 
failed  to  compute  the  leeward  side  of  £$  but  we  may  firmly  state  that  It  Is 
uncomputable  -  and  indeed  has  no  meaning  for  the  problem  originally  posed. 

Previously  several  authors,  notably  Wang  (1974)  and  Gelssler  (1976)  have 
attempted  to  compute  the  flow  pattern  in  a  region  of  moderate  or  strong  cross- 
flow  reversals  but  all  have  soon  experienced  considerable  difficulties  due  to 
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deficiencies  In  their  computational  procedure.  Wang's  calculations  are  prob¬ 
ably  the  most  successful  but  for  a  =  30°  b/a  =1/4  he  is  hardly  able  to 

) 

proceed  past  £  =  -0.88  for  e  >_  150°  the  line  of  zero  cross-flow  crossing  j 

i  at  £  ^0.92.  Using  his  knowledge  of  the  experimental  information 

about  the  laminar  boundary  layer  on  prolate  spheroid,  Wang  infers  that  the 

leeside  of  coincides  with  the  windward  side  of  *s  at  least  as  far  as 

€  =  -0.4  and  that  this  is  an  example  of  open  separation.  Our  studies 

show  that  such  a  claim  is  not  legitimate  for  it  is  quite  impossible  to 

integrate  the  equations  on  the  leeward  side  far  enough  towards  the  windward 

side  to  reach  i$.  The  concept  of  'open'  separation  has  no  place  whatsoever 

In  the  theory  of  three  dimensional  boundary  layers  on  prolate  spheroids  when 

the  pressure  gradient  Is  prescribed  once  it  Is  found  a  singularity  occurs 

at  the  ok  of  accessibility. 

There  are  cogent  reasons  for  believing  that  open  separation  Is  an 
important  feature  of  experiments  on  flow  past  thin  spheroids  at  high 
Reynolds  number.  In  addition  to  that  adduced  by  Wang,  further  evidence 
may  be  seen  In  the  beautiful  experiments  recently  reported  by  Han  and  Patel 
(1979).  However  In  all  these  examples  the  external  velocity  field  Is  not 
determined  by  Invlscld  arguments  only  but  also  involves  a  subtle  interplay  with 
the  boundary  layer.  An  Immediate  consequence  Is  that  the  singularity 
at  ts  on  the  windward  side  Is  prevented  arid  from  that  It  follows  that  t$ 
cannot  be  an  envelope  of  limiting  streamlines.  Thus  for  a  real  flow 
Llghthlll's  concept  (1963)  of  *s  as  a  limiting  streamline  Is  relevant. 

Integration  of  the  equations  across  may  well  be  possible  and  If  It  Is 
an  example  of  open  separation  we  should  be  able  to  continue  all  the  way  to 
i  (e  *  180°).  On  the  other  hand  If  Is  closed,  a  good  approximate 
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result  might  well  be  obtained  if  a  variation  of  the  Reyhner-Fliigge-Lotz 
(1968)  approach  were  used.  If  a  more  accurate  answer  were  required,  there 
are  methods  in  two  dimensions  available  for  making  use  of  downstream  condi¬ 
tions  (Williams,  1975,  Cebeci ,  Keller  and  Williams,  1979)  which  seem  capable 
of  generalization. 
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